On the turbulent sources of the solar dynamo 
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■ Abstract 

■ We revisit the possible turbulent sources of the solar dynamo. Studying axisym- 
^ , metric mean-field dynamo models, we find that the large-scale poloidal magnetic 

^1 field could be generated not only by the famous a effect, but also by the x J and 

O ■ shear-current effects. The inclusion of these additional turbulent sources alleviates 

-4— > , several of the known problems of solar mean-field dynamo models. 

^ : 1 Introduction 

> : 

^ ■ Most solar dynamo models use the scenario proposed by Parker jsj. In this scenario, 

. the solar magnetic field is produced by an interplay between differential rotation and the 

^ \ \ collective action of cyclonic turbulent convection fiows. The latter is widely known as 

^ \ the a effect which is believed to be responsible for the generation of the poloidal 

' component of the large-scale magnetic field (LSMF) of the Sun. There are, however, a 

O ■ number of problems with the a effect (H, and recent developments of mean-field magne- 

> . tohydrodynamics have stimulated a quest for additional or alternative turbulent dynamo 

^ , mechanisms. The infiuence of turbulence on the evolution of the LSMF is expressed by 

^ \ the mean electromotive force (MEMF) S = (u x b), where u and b are the fiuctuating 

■ - - ' parts of the velocity and magnetic field (angular brackets denote ensemble averages). It 
can be written in compact form as 



(d + 7) o B - (^r/ + 1] + + . . . j o (V X B) 



Here B is the mean, large-scale magnetic field. The tensors d, 7 and 17 describe the 
OL effect, turbulent pumping and turbulent diffusion, the x J effect (or b'^^^ effect) 
and W the shear-current effect (or effect). Discussions of the different contributions 
to S may, for instance, be found in @, 0|- In this paper, we use calculations of them 
as described in and construct a set of axisymmetric kinematic dynamo models in 
a convective spherical shell. These models are obtained as solutions of the mean-field 
dynamo equation 

— = V X (V X B + £) , (2) 
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Figure 1: Basic model quantities. Left panel: Contours of the rotation rate in the solar 
convection zone (left) and geometry of the pumping velocity of the toroidal LSMF (right); the 
maximum amplitude of the pumping velocity is about Im/s. Right panel: Radial profiles of the 
Coriolis number O* (solid line) and a^j)^ sec 6 (dashed line), a^^ changes sign near the bottom of 
convection zone. Similar radial dependences are found for the other components of the a tensor. 



where V = r sinOQ {r,9), with fl{r,9) modeling the differential angular rotation rate of 
the Sun (r and 9 denote radius and colatitude); see Fig. [H left-hand side in left panel. 

/A ir 9) 

The LSMF is modeled in the form B = e^B (r, 6*) + V x ^ \, (e^ denotes the 

\ rsmt^ / 

azimuthal unit vector), and the MEMF is calculated according to Eq. ([T]). Details of the 
procedure are given in [5|. The radial profiles of characteristic quantities of the turbulence, 
such as the rms convective velocity Uc, the correlation length and time ic and r^, and the 
density stratification parameter G = Vlogp (p is the mass density), are computed on the 
basis of a standard model of the solar interior j^. Some of the generation effects depend 
on the rms intensity of a fiuctuating background magnetic field h^^\ generated by a small- 
scale dynamo. We assume energy equipartition between b*^°^ and the turbulent velocity 

field, i.e. ^ (^bW^^/ (uc-\/4vrp) = 1. The relative strengths of the turbulence effects are 

controlled via parameters Ca (a effect), C^j (f2 x J effect) and Cw (shear-current effect); 
in addition, turbulent pumping is always present, depending merely on the turbulence 
level (with which all other contributions are also scaled). The integration domain is 

radially bounded by r = O.Tli?© and r = 0.96/2©, where the boundary conditions are 

drB 

— — = 0, A = at the bottom boundary, and vacuum conditions at the top boundary. 

or 

Four basic model quantities, namely, the differential rotation, the turbulent pumping 
velocity of the toroidal component of the LSMF and the radial profiles of the Coriolis 
number Q* = 2QoTc {Qq = 2.86 ■ 10~^s~^ is the surface rotation rate) and of a<^<^sec 9 
measures the strength of the azimuthal a effect) are shown in Fig. [H 



2 Results 

We now discuss examples of dynamos obtained by combining the a effect with the 17 x J 
and shear-current effects. It was found earlier [9| that for solar conditions and without 
meridional circulation, S^^^il dynamos (fi x J effect plus differential rotation) have steady 
non-oscillatory dynamo eigenmodes. We found that the same applies to S^^^^l dynamos 
(shear-current effect plus differential rotation) and to combinations of S^^^il and 5^^^^l 
dynamos. This means that the a effect is needed in solar dynamo models, at least in 
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Figure 2: Stability diagrams for the a'^d^^^Q (left panel) and a^J^^^Q (right panel) dynamos. 
The stable regions (no dynamo) lie below and to the left of the solid lines. Shading indicates 
dominance of dipole modes. 

such without meridional circulation. Fig. [2] gives stability diagrams for and 
Q,2^{iy)Q dynamos, which include the usual a'^il dynamo as a limiting case. As is seen 
in Fig. [21 in this latter case the dipole modes are not the most unstable modes. Fig. [3] 
shows snapshots of the evolution of the LSMF in the a'^fl dynamo and the associated 
butterfly diagram, obtained from the flrst unstable dipole eigenmode. We assume that the 
sunspot activity is produced by the toroidal LSMF in the whole convection zone and have 
integrated the toroidal fleld over radius. The obtained dynamo has some similarity with 
the observed evolution of the large-scale magnetic activity of the Sun, e.g., the "correct" 
phase relation between the toroidal and poloidal components — the polar reversal of the 
radial component of the LSMF takes place when the amplitude of the toroidal component 
passes its maximum, and the toroidal fields drift towards the equator in the course of 
the activity cycle. There are a number of problems, however. Two of them are: (i) The 
period of the activity cycle is T < 1 yr (corresponding to about 0.1 turbulent diffusion 
times), which is much shorter than the observed period of 22 yr. (ii) The activity maxima 
are found at mid latitudes, rather than, as observed, at low latitudes. 

The inclusion of the two additional turbulent generation mechanisms of the poloidal 
LSMF reduces the frequency of the dynamo wave. Fig. [4] demonstrates this for the example 
of an a'^S^^^^l dynamo. Oscillatory modes are excited if Ca is not very small compared 
to C^. In the range 1 < C^/Ca < 3 the match with the solar observations is best. As 
can be seen in Fig. [2l the dipole modes are the primary dynamo modes in this parameter 
interval. While the wings of the obtained dynamo waves are too wide, other qualitative 
properties, such as the drift directions of the toroidal and poloidal components of the 
LSMF and their phase relation, are reproduced correctly. The period of the dynamo in 
this parameter range is shorter than but comparable with the turbulent diffusion time. 

Our last example is a dynamo where the poloidal LSMF is built up with the help of 
the shear-current effect. Fig. [5] shows the time evolution of the LSMF and the associated 
butterfly diagram for an a'^S^^^^l dynamo model. Some properties of the model agree well 
with the observations. For example, the period of the dynamo waves is comparable with 
the turbulent diffusion time. Also, the maximum amplitudes of the toroidal component 
of the LSMF are found below 30° latitude, and this fleld component drifts towards the 
equator during the cycle. A very unwelcome issue of the model, however, is that the 
maxima of the radial component of the LSMF are also concentrated near the equator. 
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Figure 3: a^Q model on the basis of the first unstable dipole eigenmode. Top: Snapshots of 
the strength of the toroidal LSMF (greyscale plot) and field lines of the poloidal LSMF over 
an half cycle. Bottom: Butterfly diagram in the form of contour lines of the toroidal LSMF 
(integrated over radius) with overlaid greyscale plot for the radial LSMF at the top boundary. 
Time is measured in units of the turbulent diffusion time. 
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Figure 4: As Figure El but for the a'^d^^'^n model with = 0.1, Ca = 0.05. 
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Figure 5: As Figure O but for the a'^d'-^^n model with Cw = 0.1, C„ = 0.01. 

3 Conclusions 

The inclusion of the fix J and shear-current effects in addition to the a effect and differ- 
ential rotation in axisymmetric kinematic dynamo models alleviates some of the known 
problems of current mean-field solar dynamo models. For instance, the simulated period 
of the dynamo is increased and comes, thus, closer to the observed activity period. Fur- 
thermore, the large-scale toroidal field is concentrated towards the equator, thus bringing 
the models in better agreement with the observation of two activity belts relatively close 
to the equator. Improved diagnostic tools based on both observation and theory are 
needed to clarify the roles of the different turbulence effects, as well as that of meridional 
circulation, which was not included in the present study. 
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